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Abstract 



^ . We calculate the 0{aas) corrections to the relationships between the MS Yukawa 

iy~, I couplings and the pole masses of the first five quark flavours in the standard model. 

CNl ' We also present the corresponding relationships between the MS and pole masses, 

^^ . which emerge as by-products of our main analysis. The occurring self-energies are 

.^ I evaluated using the method of asymptotic expansion. 
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1 Introduction 

One of the most intriguing puzzles in contemporary elementary particle physics is related 
to the deciphering of the seemingly arbitrary pattern of the fermion masses and the 
elements of the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix of the standard 
model (SM) in the framework of some more fundamental theory. It is generally believed 
that grand unified theories (GUTs), possibly realized in the context of supersymmetry, are 
able to provide a key to the understanding of this fundamental problem |lj. The crucial 
idea is that a judiciously chosen set of independent parameters, appropriately evolved to 
the GUT scale, obey simple relationships. Essential ingredients of this formalism include 
renormalization group equations, which determine the scale dependence of the running 
parameters, and threshold relations, which relate the running parameters at the scales of 
the elementary-particle physics experiments to the physical parameters, e.g. pole masses 
and CKM matrix elements, extracted from the latter. It has become customary to define 
the running parameters in the modified minimal-subtraction MS scheme. 

The Callan-Symanzik beta functions of the SM and its most popular extensions are 
well established at one loop and beyond [2]. As for the threshold relations of the Yukawa 
couplings and CKM matrix elements the current status is as follows. The relationships 
between the MS and on-shell definitions of the CKM matrix elements are known at one 
loop ^. The relationships between the MS and pole definitions of the quark masses 
have been elaborated at one 4J, two j3], and three jH] loops in quantum chromodynamics 
(QCD). These pure QCD corrections readily carry over to the relationships between the 
MS Yukawa couplings and the pole masses of the quarks, which are actually relevant for 
GUT analyses. However, the situation becomes more involved if electroweak quantum 
corrections are taken into account. Then, the relationships between the MS definitions of 
the fermion masses and the Yukawa couplings become nontrivial [HIHI. The full one- loop 
corrections to the relationships between the MS Yukawa couplings and pole masses of the 
SM fermions were derived in Ref. [f . They were found to be gauge-parameter independent 
and devoid of tadpole contributions. However, in that reference, also a gauge-parameter- 
independent MS definition of fermion mass was presented, and tadpole contributions were 
found to be essential ingredients for that. 

In this paper, we take the next step and evaluate the mixed two-loop corrections, 
involving one power of Sommerfeld's fine-structure constant a and one power of the strong- 
coupling constant ««, to the relationships between the MS Yukawa couplings and the pole 
masses of the first five quark flavours in the SM. For simplicity, the calculation is performed 
in the limit where the third quark generation does not mix with the first two, i.e. where 
Vub = Vcb = Vid = Vts = 0, which is approximately realized in nature p. As in Ref. jTj, 
we pay special attention to the tadpole contributions. Again, they are found to cancel 
in the relationships between the MS Yukawa couplings and the pole masses, while they 
are indispensable to render the MS definitions of quark mass independent of the gauge 
parameters. 

As for the top quark, the 0(aas) correction to the relationship between the MS and 
pole masses was recently evaluated in Ref. JU] retaining the full mass dependence. Also 



there, the tadpole contribution was found to be necessary in order to get rid of the 
gauge-parameter dependence. An alternative definition of MS top-quark mass, without 
inclusion of the tadpole contribution, was proposed in Ref. ^l]. Since Ref. ^Ij is dealing 
with leading heavy-top-quark effects, the gauge-parameter dependence does not yet show 
up. However, we caution the reader that it will once subleading corrections are to be 
included in this formalism. Then, the formulas that express physical observables in terms 
of the so-defined MS top-quark mass will explicitly depend on gauge parameters as will 
the value of this mass. Although this is not forbidden by first principles, it is cumbersome 
in practice because, whenever a value of this mass is to be quoted, the underlying choice 
of gauge needs to be specified as well [7|. 

This paper is organized as follows. In Section |2l we set up the theoretical framework 
and derive a master formula for the relationship between the MS Yukawa coupling and 
the pole mass of a quark through 0{aas)- In Sectional we briefly recall the method 
of asymptotic expansion and apply it to recover the 0{a) results for the first five quark 
flavours. We then illustrate the goodness of the resulting expansion for the case of bottom 
by comparison with the analytic 0{a) result of Ref. [7j. The 0{aas) results for the first 
five quark flavours are presented in Section HJ As a by-product of our analysis, we also 
find the relationships between the MS and pole masses of the first five quark flavours to 
this order. They can be found in Section |S| Finally, we conclude with a summary in 
Sectional 

2 Deriving the relations 

The relationship between the MS Yukawa coupling /i/(/i) at renormalization scale /i and 
the pole mass Mf of a fermion / can generically be written as 

hj{f^)=2'/'GTMf[l + 6j{f,)], (1) 

where Gp is Fermi's constant. Here, we compute the corrections Sf(^) for the first five 
quark flavours to 0{a), 0{as), and 0{aas)- 

In the SM, the fermion masses are generated from the Yukawa interactions through 
the Higgs mechanism of spontaneous electroweak symmetry breaking. In terms of bare 
parameters, which are henceforth labelled by the subscript zero, this is manifested through 
the identity 

T^ffi = "I^Z-o. (2) 

where v is the vacuum expectation value of the Higgs field. Next, we relate the bare 
quantities in Eq. Q with their renormalized counterparts in Eq. (^. 

The pole mass is defined as the zero of the inverse propagator. The inverse fermion 
propagator can be written as 

iSfiqy^ =^ -nif^Q + J^fiq), (3) 
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Figure 1: Sample diagrams which contribute to the bottom-quark self-energy. 
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is the self-energy of fermion /. In Eq. (JH), the subscripts V, A, and S denote the vector, 
axial- vector, and scalar components, respectively, and the sum runs over all contributions 
through the desired order. Here, we consider the one-loop contributions of 0{a) and 
0{as) and the two-loop ones of 0{aas), which we denote by the superscripts (1) and 
(2), respectively. Some sample Feynman diagrams are depicted in Fig.[Tl^a)-(c). Inserting 
Eq. (JH) into Eq. (jH)) and setting the latter to zero leads to jH] 
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{l - S«(l) - Sf (1) + 4^)(1) [S« (1) - 24'^'(1)] } 
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which is valid through the two-loop order. Here, the prime indicates differentiation with 
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respect to rrifQ/q, and the notation SV = Tjty + ^j s has been introduced. Note that, 
while deriving Eq. (J3J, an expansion of tti/q about Mf has been performed in order to 
obtain unity as the argument in the self-energies. Furthermore, we remark that EI^ does 
not appear in Eq. (j3)). For later convenience, we also provide the inverted relation, which 
reads 
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^(1) 
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By construction, the mass of fermion / which appears in the self-energies on the r.h.s. of 
Eqs. (0) and © is the pole mass. These equations become manifestly gauge-parameter 
independent if one includes the tadpole contributions Tj in the self-energies j7], i.e. if the 
tadpole- free self-energies S^ in Eqs. (0) and (0) are replaced by S^ +TJ . Some sample 
tadpole diagrams are shown in Fig. ^d)-(f). 



The Fermi constant Gp is not a basic parameter of the SM Lagrangian density. This 
means that it has to be expressed in terms of such parameters. To do this, one calculates 
a process in the SM and in the Fermi theory, and equates the two results. This was first 
done in a classical paper by Sirlin |^, who considered the muon lifetime. Adapting his 
result for corrections of 0{aas), we have 
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Here, awwil"^) is the tadpole- free W^-boson self-energy at four- momentum q, and Tww 
denotes the corresponding tadpole contribution. The quantity E contains those wave- 
function renormalization, vertex and box corrections to the muon decay width which the 
SM introduces on top of the Fermi model improved by quantum electrodynamics (QED). 
At 0{a) and in 't Hooft-Feynman gauge, it may be written as 
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where we have used the abbreviation c?„ 



M^/Ml 



Here and in the following. 



we use dimensional regularization, with d = 4 — 2e space-time dimensions and 't Hooft 
mass scale fi, adopt a MS-like notation where the typical combination 'Je ~ ln(47r) is 
suppressed, and do not display terms of 0{e). There are no corrections of 0{aas) to this 
quantity. At the one-loop level, the M^-boson self-energy can be split into a bosonic and 
a fermionic part as 

Uwwiq') = Kwil") + ^wwiQ"), (9) 

and, in the 't Hooft-Feynman gauge, we have (see, e.g., Ref. [7j) 
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where A^c = 3 denotes the number of quark colours and we have omitted terms quartic in 
rrib^o in Eq. |TT|. Note that only H^p^/(0) receives 0{aas) corrections. 

Finally, the relation between the bare and the MS-renormahzed Yukawa couplings is 
given by 
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where Shj, Sj^qj^, and (5}^cr denote the appropriate MS counterterms. 



Inserting Eqs. ©, ©, and (fT^ in Eq. Q and comparing the result with Eq. (P), 
where we decompose Sf{fi) = 6j (yu) + 6j (/i), we obtain the master equations 
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Here, all self-energies are tadpole-free and it is understood that all terms of higher orders 
in the electroweak couplings are discarded. We emphasize again that Eqs. (jHI) and ((Tj) have 
to include tadpole contributions to ensure their gauge-parameter independence. However, 
in Eqs. fjl3|) and ()14|) . these contributions cancel thanks to the identities 
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which relate the tadpole contributions to the fermion and boson self-energies at 0{a) 
and 0{aas), respectively. We have verified Eq. (fT^ for arbitrary values of the gauge 
parameters. Note that Eqs. (fr5|) and (fT^ are finite and gauge-parameter independent. 
The later follows immediately from the gauge-parameter independence of Eqs. (jH)) and 
((7j). It should also be remarked that the 0(e) terms of the quantities Hvi/h/(0) and E are 
not needed explicitly, as they drop out in the combination contained in Eq. (fT^ . 
Through 0{aas), 5f[ii) may be decomposed as 
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where Cp = {N^ — l)/{2Nc) is the eigenvalue of the Casimir operator of the fundamental 
representation of SU(3)c. The well-known one-loop QCD contribution reads: 
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The two-loop [5] and three-loop jH] QCD contributions are also known. The 0{a) con- 
tribution is given in Ref. [7j. The evaluation of the 0{aas) one is the main purpose of 
this paper. It is outlined in Section HI Due to the presence of many different mass scales, 
an exact calculation is rather cumbersome. Thus, we adopt the concept of asymptotic 
expansion, which we recall and apply to the one-loop case in the next section. 



3 Asymptotic expansion 

The method of asymptotic expansion y^ allows for the evaluation of Feynman integrals 
which depend on various parameters of different sizes. This is achieved by expanding the 
integrand in small ratios of these parameters. As a result, one obtains an infinite series 
in powers and logarithms of the expansion parameters. 

In our case, the so-called hard-mass procedure is used, where one mass M is much 
larger than all other masses {rrii} and external momenta {qi}- In this limit, the expansion 
can be performed in a diagrammatic way, according to the prescription 

r(M, m, q) ''^°° Y. r/7(g, m) * 7^g„™,)7(M, m^, g^), (18) 
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where F is the initial diagram and the operator 7(q^^m^) performs a Taylor expansion in all 
small parameters of the subdiagram 7. The sum runs over all subdiagrams that contain 
the large mass M and are one-particle-irreducible in all connected parts after contracting 
all lines with the large mass. 

Let us first consider the case of bottom. In terms of the inherent mass scales, the 
bottom-quark self-energy diagrams can be divided into two classes: 

1. those involving the exchange of charged bosons (PF, 0), and 

2. those involving the exchange of neutral bosons (Z, %, H). 

Exploiting the fact that M^ = ^/^Mw and M^ = \f^Mz depend on the arbitrary gauge 
parameters ^vk and ^z-, we can arrange for the masses involved in these diagrams to fulfil 
the following hierarchies: 

Mfe < Mw, M^ < Mt, (19) 

Mh<t:Mz,M^,MH. (20) 

Since our final results are gauge-parameter independent, we can replace Eq. (jT^ by the 
alternative hierarchy^ 

Mb<^Mw <^Mt<^ M^. (21) 

This changes the intermediate results for some of the diagrams, but not the final result 
for 6b{fJ.). 

In the case of the light quarks q = u,d, s, c, the situation is much simpler. Since we 
neglect their masses nig, except for the linear appearances in front of ^q,s, all diagrams 
are either reduced to vacuum integrals where the scale is given by a boson mass or to 
on-shell integrals which can be adopted from the bottom-quark case. 

Similarly, in the case of the fermionic part of the PF-boson self-energy, which has to 
be evaluated for vanishing external momentum, we only have the hierarchy Mh <^ Mt. 

^In principle, we could also choose M^,M^ <C Mj,. In this case, however, Eq. (|18|l could not be 
applied, since the external momentum is on the bottom-quark mass shell, which requires different rules 
of calculation J^! ■ 



As for the bosonic part of the W^-boson self-energy, we use the exact expression given in 

Eq. (uni). 

The calculation of the self-energy diagrams can be performed in a completely auto- 
mated way. This is achieved by the successive use of the computer programs QGRAF [1^, 
q2e ^, exp [IE], and MATAD ^ (for a review, see Ref. |IH|). First, QGRAF is used to 
generate the Feynman diagrams. Its output is then rewritten by q2e to be understandable 
by exp. This program performs an asymptotic expansion of the diagrams, if necessary, 
in an iterated fashion. Finally, the FORM-based |19j program MATAD is used to perform 
the actual calculations and expansions in e. For the calculation of certain diagrams, it is 
supplemented with the package 0N-SHELL2 pOj . 

In a first step, we explicitly check the gauge-parameter independence of the combina- 
tion Sj (1) + Tr appearing in Eq. (jH)). We refrain from listing the corresponding results, 
which can be extracted from the expressions given in Section where the relation between 
the MS and pole masses is discussed. 

Let us now move on to 5^° . Since E is only known in 't Hooft-Feynman gauge, we 
have to choose this gauge also for S^ (1) and nyy(/(0). We determine the counterterm 
in Eq. (fT^ by requiring that Eq. (fTSj) be finite and thus find 



r(i) 

Of,CT 



^b,CT 



as 3 1 a 

OiT- - + - 

n 4 e 7r 



,L;2 _ 3^2 _ ]_ 

4 4 s^ 



+ 



771 



t,0 



M'w 



1 
I2" 



3 

-32 + 



16 




M'w 



(22) 



Inserting all the above ingredients into Eq. ()13|). we obtain 
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Here and in the following, ellipses stand for terms of 0{M^/M^), 0{M^,M^/M^), or 
0{M^ /M^). As we will see below, they are very small and can safely be neglected. Note 



that there are no terms proportional to M'^ /M^ in Eq. ()23p. Here and in the following, Qf 
is the fractional electric charge of fermion /, I-^j is the third component of weak isospin of 



its left-handed component, Vf 



{lf-2slQf)/{2c^s, 
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If /{2c, 



•Sw) are its vector 



and axial-vector couplings to the Z boson, respectively, and we use the abbreviations 
h = Hf^VMf) and hj = ln(M2/M|), 

For the light quarks q = u,d,s,c, we have 
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Let us now discuss the quality of our approximation by comparing 5]^ (/x) of Eq. PH|) 
with the exact result obtained in Ref. |7|. As input values for our numerical analy- 
sis, we use a = 1/137.035, Mb = 4.5 GeV, Mt = 174.3 GeV, Mw = 80.423 GeV, 
Mz = 91.1876 GeV 0, and Mh = 120 GeV. Furthermore, we choose /z = Mb. Fig- 
ure |21 displays our result as well as the exact one as functions of M^/Mt. To illustrate 
the convergence of the expansion, we show the leading-order contribution separately and 
then successively add the subleading terms. It is obvious from this figure that the se- 
ries in M^r/M^ converges very fast. Indeed, in the range Mw/Mt ^0.8, the sum of all 
calculated terms can barely be distinguished from the exact result. The approximation 
including the 0{M^ /M}^) terms actually provides an excellent approximation also for 
Mw/Mf ^ 1.5, whereas the lower-order approximations start to exhibit significant devi- 
ations at Mw/Mt ^ 1. 

4 Two-loop result and numerical analysis 



In order to obtain (5/(/i) to 0{aas), one has to evaluate the respective two-loop corrections 
to S/(g) and ni^^(O), as can be seen from Eq. (fT^ . The result for the latter is well-known 
(see, e.g., Ref. |2I1)- Including terms quadratic in m?,, one has 
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where C denotes Riemann's zeta function, with the value ((2) = tt^/6. We again refrain 
from listing explicitly our results for the fermion self-energies, but provide the results 
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Figure 2: Comparison of our result for (a/'Tr)^^" (M;,) with the exact result. The latter is 
shown as a continuous line. The long-dashed line indicates the leading-order contribution 
proportional to M^ and M^ only. The contributions given by successively adding the 
subleading terms are indicated by the short-dashed lines. The horizontal line marks the 
actual value of Mw/Mf. Note that for this plot Mw is fixed to its default value and Mt 
is varied. 



for Sr (jj,), from which the expressions for Hr \ 

v(2) 



1) can be obtained if the result for the 



counterterm 6f Ijj^ introduced in Eq. (fT2j) is given. Similarly to the one-loop case discussed 



in the previous section, Sj^ct i^ determined by requiring that Eq. (jT^ be finite upon 
renormalization of the top- and bottom-quark masses, and we find 
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Note that there are no poles proportional to mff^ or m'^Q, although such terms do appear in 
certain tadpole diagrams. However, the tadpole contributions cancel in the computation 
of Sb{n), as we have explained in Section |2l As an additional check, we can exploit the 
fact that the non-local logarithmic terms have to cancel as well. 

Finally, the 0{aas) correction to the relationship between the MS Yukawa coupling 
and the pole mass of the bottom quark turns out to be 
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We should mention that the QED portion of the 0{aas) contribution can also be inferred 
from the two-loop 0{a1) one j3] by replacing a'lC'jp with 2asCF(y.Q'^f and setting all the 

other colour structures to zero. In order to obtain the result for 5^°" {jj) in Eq. (f?fjl . we 
have used the tadpole-free bottom-quark self-energy and adopted the 't Hooft-Feynman 
gauge, the one in which E is given. The gauge-parameter independence of the full self- 
energy is discussed in the next section. 

In the case of the light quarks q = u,d, s, c, we obtain the compact expression 
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In order to demonstrate the good convergence properties of our result, we show in Fig.|Hl 
the 0{aas) contribution to 6b{Mf)) as a function oi Mw/Mt for Mh = 120 GeV. In addition 
to the parameters specified after Eq. (j^ . we use as{Mb) = 0.1905 appropriate for six 
active quark flavours, which we evaluate from the present world average a^^\Mz) = 0.1172 
[H] using the program RunDec |23- Again, we observe an excellent convergence up to 
Mw/Mt ~ 0.8, which is more than sufficient for our purpose. Since it is not possible to 
distinguish between the different lines in the physically interesting region in Fig. IHl we 
give explicit numbers for the subleading terms at Mw/Mt ~ 0.46. For this purpose, we 
write the series as 
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Figure 3: (aa./Tr^)^^^""'^!^^) as a function of Mw/Mt for Mh = 120 GeV. The 
long-dashed hne indicates the leading-order contribution only. The contributions given 
by successively adding the subleading terms are indicated by the short-dashed lines. The 
horizontal line marks the actual value of Mw/Mt. Note that for this plot Mw is fixed to 
its default value and Mt is varied. 



and thus obtain 



2\m 



m=0 

= 0.997925 + 0.004058 - 0.001451 - 0.000448 - 0.000078 - 0.000006. 



(30) 



This demonstrates that, aside from the leading-order terms, the (^(M^/Mj') term yields 
the largest contribution by far, while all other subleading terms are rather small. Fur- 
thermore, we wish to mention that the inclusion of the corrections quadratic in M^ have 
no visible effect in Fig. El 

Our final numerical results are presented in Figs.EJandEJ Figure|^shows the 0{a) and 
0{aas) contributions to Sb{Mb) as functions of M^. For comparison, also the pure QCD 
ones of 0{as), C(a^); and 0{a^) are plotted. We observe that the 0{aas) contribution 
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Figure 4: 6b{Mh) as a function of Mh- The dash-dotted and full lines indicate the 0{a) 
and 0{aas) contributions, respectively. For comparison, also the pure QCD contributions 
of 0(q;s) (short-dashed line), O^a"^) (long-dashed line), and 0{a^) (dotted line) are shown. 
The latter are, of course, independent of Mh- 



exhibits a rather weak dependence on Mh- Futhermore, for Mh ~ 500 GeV, it becomes 
comparable in size to the 0{a) contribution, whose magnitude decreases with increasing 
value of Mh- From Fig. HI we also see that our new 0{aas) contribution is of the same 
order of magnitude as the 0{al) one. 

Figure El shows ^^(/i) for Mh = 120 GeV as a function of /i in the range M^ < fi < 
1000 GeV. We observe that the 0{aas) contribution takes positive values in the range 
20<,fi<,600 GeV. It is comparable in size to the 0{a^) one only in the lower /i range. 



5 Relationship between MS and pole masses 

As a by-product of our calculation, we obtain the 0{aas) correction to the relationship 
between the MS and pole definitions of mass for the first five quark flavours. As dis- 
cussed in Section 121 it is necessary to include the tadpole contributions in order to obtain 
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Figure 5: 6b{fi) as a function of /i for Mh = 120 GeV. The same coding as in Fig. HI is 
adopted. 



gauge-parameter-independent results. This has been checked exphcitly by considering the 
hierarchies discussed in Eqs. (IT^ . (|^. and (PT|) . 

The starting point for the derivation of the MS to on-shell relationship is Eq. 0, 
where we replace the bare quark mass by the MS one. The corresponding relation is 
obtained in analogy to Eq. (j^Uj) . with the only difference that we also have to allow for 
terms quartic in rrib^Q and mtfi because of the tadpole contribution. By requiring finiteness 
of the resulting relation, we obtain 



nib 



1 



ry. 




3 1 




cy 




Cf— 


+ 




TV 




4e 




TT 



M'w 



84 M^ 






--AQI + 

N, 



v^-a^]- 



2 



324 M^ 



+ 



rn; 



Ml 



IQclsl M^ 



9 
16 



Q 



44 Ml,Ml 
9 



+ 



ml 



324 M2, 



Ml 



Ml 



H 



324 M^ 



Ml 



6444 Ml 



1284 Mir 



15 



9 

324 


Ml 


+ 


9 m? 

644 Ma^ 


3 
32 


{Ql 


+ 


^t) 


21 2 

+ 32^^^ 


3 

324 




+ 


84 


ml 
M^Ml 



9Nc fn^ 9 ml 

164 Mp^ ~ 644 M|: 
9 



1284 
3 



+ 



ml 



324 M^ 



+ 



(31) 



If the tadpole contributions were omitted here, then there would be no terms quartic 
in rrib and fnt and no terms involving Mh- In turn, gauge-parameter dependence would 
appear in the nonleading terms. 

The relationship between the MS and pole masses of the bottom quark through 0{aas) 
reads: 
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At 0{aas), the term quadratic in Mt agrees with the result of Ref. 

We conclude this section by presenting the relationship through 0{aas) between the 
MS and pole masses of the first four quark flavours. It reads: 
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6 Summary and conclusion 

We calculated the 0{aas) corrections to the relationships between the MS Yukawa cou- 
plings and the pole masses of the first five quark flavours in the SM. We demonstrated that 
these relationships are gauge-parameter independent and free of tadpole contributions. 

The method of asymptotic expansion was used in the calculation, and the results 
were found as polynomials in small mass ratios with coefficients that contain logarithms 
of these mass ratios. The goodness of this method was tested in two ways. On the one 
hand, our 0{a) expression for 6b{Mb) was found to reproduce the analytic result of Ref. [7| 
extremely well. On the other hand, our 0{aas) result turned out to converge very fast 
as higher powers of M^/M^ were included. These two observations reassure us of the 
reliability of our result, which is equivalent to the exact analytic result for all practical 
purposes. 

As for the phenomenological significance of our results, our numerical analysis revealed 
that, at Mh = 100 GeV, the 0{aas) contribution to Sb{Mb) amounts to roughly one half 
of the 0{a) one. These two contributions depend very differently on the value of Mh, 
and they even cross over at Mh ~ 500 GeV. Obviously, the 0{aas) contribution cannot 
be neglected against the 0{a) one. 

We also calculated the 0{aas) corrections to the relationships between the MS and 
pole masses of the first five quark flavours. We showed that these relationships are gauge- 
parameter independent if the tadpole contributions are properly included in the contribut- 
ing self-energies. 
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